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We show that applying strain on half-doped manganites makes it possible to tune the system 
to the proximity of a metal-insulator transition and thereby generate a colossal magnetoresistance 
(CMR) response. From our realistic microscopic model of strain and magneto-transport calculations 
within the Kubo formalism, we demonstrate that under tensile strain a ferromagnetic charge ordered 
insulator, previously inaccessible to experiments, becomes stable. The phase competition not only 
allows control of CMR in ferromagnetic metallic manganites but further generates CMR response 
in otherwise robust insulators at half doping. 



Introduction. — Transition metal oxides have long been 
studied for their surprising emergent behavior such as 
high T c superconductivity in the cuprates, ferroelectricity 
in the titanates, and colossal magnetoresistance (CMR) 
in the manganites. However, very recent advances in het- 
eroj unction growth [1-3] have opened the possibility of 
producing atomically perfect interfaces of oxide materi- 
als, and therefore applying precisely controlled strain and 
otherwise tuning the properties of oxide thin films. In 
this Letter we address the question of how strain impacts 
CMR in the manganites, in particular at "half-doping." 
The materials in this class have a prototypical chemical 
formula of Ao.sA'o.sMnOa where A is a rare earth and A' 
an alkaline earth metal [4]. At large bandwidths (BW), 
the half-doped manganites are ferromagnetic metals (F- 
M), while narrow BW materials are spin, charge, and 
orbitally ordered insulators (SCO-I) as described later in 
the text. CMR is known to occur in F-M materials close 
to the metal-insulator phase boundary, and is the result 
of phase competition (between the F-M and a charge- 
ordered insulator), which is traditionally controlled by 
isovalent chemical substitutions at the A-site. [4] 

Most prior work on the effect of strain (both theory [5- 
11] and experiment [12-20]) has focused on how it affects 
the magnetic and electronic phases, with little emphasis 
on magnetotransport. In this Letter, we propose a real- 
istic microscopic model of the manganites under strain 
and use it to calculate magnetoresistance at finite tem- 
perature within the Kubo formalism. In addition to de- 
veloping a broader theoretical understanding of strain in 
correlated systems, we obtain the following results: 

(i) Tensile strain can stabilize a ferromagnetic charge 
ordered insulator (FC-I). This phase has not been con- 
clusively observed in any half doped manganite with tol- 
erance factor variations, but should be observable with 
strain engineering. 

(ii) We demonstrate that strain can induce phase tran- 
sitions. As a consequence, we show that the CMR in 
F-M materials can be enhanced by tuning the proxim- 
ity to the metal-insulator transition, and that insulating 
phases can be made metallic under strain and therefore 



also exhibit a CMR response. This greatly expands the 
family of materials with potential device applications. 

Model. — We begin with the "standard model" for the 
manganites. Because of the octahedral crystal field, the 
Mn e g levels have a higher energy than the t<i g levels. 
Combined with a large Hund's coupling, which ensures 
that the electron spins align ferromagnetically, this local- 
izes three Mn 3d electrons in the ti g levels which form 
local S = 3/2 moments (called "core spins"). The re- 
maining electrons, if any, are itinerant and occupy two 
bands that result from the hybridization of the Mn e g 
levels. The model also includes an effective antiferro- 
magnetic superexchange J between neighboring Mn core 
spins, and finally the e g electrons couple to Jahn- Teller 
phonons with a coupling strength A. 

All energy scales are given in units of the unstrained 
bandwidth t ~ 0.3eV. This Hamiltonian yields a remark- 
ably accurate description of the physics of the mangan- 
ites [21] and is discussed in detail in the Supplementary 
Information section I. We now extend this model to in- 
corporate the effects of strain. 

Fig. 1(a) schematically illustrates substrate-induced, 
in-plane tensile strain. Tensile (compressive) strain is 
caused by growing the film on a substrate with lattice 
parameters larger (smaller) than those of the unstrained 
film. We assume that strain is applied parallel to the 
(a-b) plane, which we take to be the plane of the d x 2_ y 2 
orbital. We quantify strain by the parameter ey = [(a s — 
a) /a] = 5a/ a, where a s and a are the substrate and film 
lattice parameters. Here a is the distance between two 
nearest neighbor Mn atoms. 

We consider volume conserving strain and use the re- 
lation e± = — 4z/ey where v is the Poisson ratio. We 
choose v = 0.375, consistent with previous estimates that 
v is between 0.3-0.4 [22, 23] for the manganites. In-plane 
compressive strain corresponds to ey < 0, while ey > 
for tensile strain. 

With this in mind, we propose three microscopic effects 
of strain that must be included in our model: 

(i) Strain modifies the hopping matrix elements. Strain 
affects the lattice parameter a which in turn modifies 
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FIG. 1. (a) Uniform tensile strain in the a-b plane is in- 
troduced by lattice mismatch with the substrate. Volume 
conserving tensile strain expands the in-plane Mn-0 bonds 
while contracting bonds in the c direction. This leads to a 
higher occupancy of the in-plane d x i_ y i orbital, a larger out- 
of-plane hopping and smaller in-plane hopping compared to 
the unstrained values. Compressive strain, not shown, has 
the opposite effect, (b) Schematic of a generic Mn-O-Mn 
bond under strain. We indicate the Mn-0 bond length d, 
Mn-O-Mn bond angle <j) and the shortest Mn-Mn distance, 
the lattice parameter a. An expansion of a along the green 
arrows causes a change in d and (j). 

both the Mn-0 bond length d and the Mn-O-Mn bond 
angle (j) as seen in Fig. 1(b). The bandwidth depends 
on both d and (j) directly through the p — d overlap in- 
tegral. In the Supplementary Information section I (c), 
we show that the effect on the hopping matrix elements 
due to the change in <\> can be neglected for a large class 
of half doped manganites, and that the hopping in the 
(a-b) plane scales with strain as t xy —> t xy (l — Ten). We 
restrict our calculations to a single layer manganite film 
in the a-b plane (as depicted in Fig. 1(a)) and refer to 
the unstrained in-plane hopping parameter t xy as £, and 
under strain as t. 

(ii) Strain modifies the antiferromagnetic superex- 
change. The superexchange coupling also scales with the 
hopping. Using considerations similar to those above 
it is shown that the in-plane superexchange scales as 
Jxy (J\\/t)0- ~ 14e || ) with strain. We will refer to 
the unstrained in-plane superexchange as J and in the 
strained case, J. 

(iii) Strain generates an orbital bias. Because of 
the increase in the in-plane Mn-0 bond length, ten- 
sile strain makes occupation of the d x 2_ y 2 orbital en- 
ergetically favorable. In-plane compressive strain favors 
the out-of-plane d 3z 2_ r 2 orbital. This orbital bias in- 
duced by in-plane compressive and in-plane tensile strain 
in Lao.7Sro.3Mn03 has been observed in x-ray absorp- 
tion [24] as well as in angle resolved photoemission [25]. 
We incorporate this in our model Hamiltonian with an 
extra term, H hias = ^Zi, a e ^ n i,a where e a = 5/2 for 
a = d% z 2_ r 2 and e a = —5/2 for a = d x 2_ y 2. 

From experiments, the e g splitting has been estimated 



to be between OAt and 3t [24, 26]. To perform calcu- 
lations, we make a conservative estimate and take the 
magnitude of the bias to be 5 ~ lOeyt, i.e., a splitting of 
±0.2£ for ±2% strain. These values are consistent with 
density functional estimates [27, 28]. 

As mentioned, we perform our calculations in two di- 
mensions, describing a single-layer manganite film in the 
a-b plane. Further, we assume strain to be uniform in the 
layer. This is sufficient to bring out the important fea- 
tures of the phase diagram and in-plane transport. We 
describe in detail our method of solution and calcula- 
tions in Supplementary Information sections II and III. 
We focus here on our results. 

Strain driven phase transitions. — In Fig. 2(a) and (c) 
we show the effect of strain on a typical large BW F-M 
and small BW SCO-I respectively. The SCO-I is an in- 
sulator with planar checkerboard charge order (CO), al- 
ternating d x 2_ r 2 j d y 2_ r 2 orbital order (00) (on the sites 
with larger charge density) , and CE type spin order (zig- 
zag ferromagnetic chains coupled antiferromagnetically) . 
In (b) we show the T=0 A — J phase diagram. On this 
we denote two unstrained parameter points, F-M (blue 
dot) and SCO-I (red dot). The results of strain on these 
are shown in (a) and (c), respectively. 

Starting from unstrained parameters, in-plane com- 
pressive strain favors F-M as seen in both (a) and (c). 
There are two competing effects here, however. Com- 
pressive strain increases the in-plane hopping which in 
turn reduces A/f, as seen by following the dashed black 
arrow in (b). This favors a metallic state where double- 
exchange promotes ferromagnetism. On the other hand, 
compressive strain increases j/t, which tends to narrow 
the BW, while the orbital bias promotes occupancy of 
the out-of-plane d 3z 2_ r 2 orbital. Both of these latter ef- 
fects work against the stability of the F-M, but we find 
the F-M to be dominant up to the maximum values of 
strain we have considered. 

In-plane tensile strain stabilizes insulators that can 
have either ferromagnetic or antiferromagnetic spin tex- 
tures as seen in (a) and (c) respectively. The insula- 
tor favored by tensile strain has long-range checkerboard 
charge order. It is stabilized by the reduced in-plane hy- 
bridization and consequently increased A/t, as seen by 
following the grey dashed arrow in (b). This effect tends 
to localize the electrons. While sufficient increase in A 
eventually turns the system insulating, regardless of the 
unstrained F-M parameter, the magnetic order depends 
crucially on the value of J/t as seen in (b). So the final 
magnetic state depends both on the amount of strain and 
the initial unstrained starting point. 

This dependence of the magnetic/charge-ordering 
scales on strain has been seen in experiments, both 
away from [15] and at half-doping [13, 19]. They 
include increasing T c with compressive strain in 
Lao.8Bao.2Mn03, suppressing T c with tensile strain in 
Lao.67Cao.33Mn03, and increasing Too with tensile 
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FIG. 2. (a) The ferromagnetic metal (F-M), ferromagnetic charge ordered insulator (FC-I), and the paramagnetic insulator 
(P-I) in the temperature (T) vs. strain (en) plane. Compressive strain (en < 0) on F-M raises T c , while tensile strain initially 
suppresses T c and beyond ~ 1% strain the system is driven into an FC-I. (b) Unstrained T — phase diagram in the A/£, 
J/t plane. This phase diagram also shows the spin-charge-orbital ordered insulator (SCO-I) and the A- type antiferromagnetic 
metal (A-M) [29]. The blue dot at (A/£, J/t) = ( 1.6, 0.08) refers to the unstrained starting point from which we calculate (a), 
while the red dot at (A/£, J/t)=( 1.65, 0.09) refers to the unstrained point from which we calculate the phase diagram in (c). 
Under strain, the hopping parameter t is rescaled which changes both the ratios A/t and J/t. The symbols at the arrow heads 
indicate such values for 2% compressive (triangles) and tensile (squares) strain respectively, (c) Tensile strain on the SCO-I 
enhances the charge ordering temperature Tco while compressive strain causes a SCO-I to F-M transition. Tce denotes the 
spin ordering temperature. In (a) and (c), for clarity, we do not show strain-induced equilibrium phase separation at low 
temperatures as this is not central to our results here [30]. 



strain on Pro.sCao.sMnOs [13, 19]. 

Stability of the FC-I phase. — From Fig. 2 (b) we see 
that suitable amounts of tensile strain on F-M with 
J/t~0.05-0.08 would convert the system into a FC-I, just 
as that depicted for (X/t = 1.6, J/t = 0.08) by the grey 
dashed arrow. This FC-I phase, in the unstrained phase 
diagram, was predicted in theory [31] at A, J values as 
in Fig. 2 (b). In small BW half-doped manganites such 
as Lao.sCao.sMnOs, signatures of this phase coexisting 
with AF-CO phase were reported [32] at 90K. This im- 
plies that in the half-doped manganites either FC-I is 
the true ground state only in a narrow range of the A — J 
phase diagram or it is a metastable state with energy 
very close to the true ground state. We predict that ten- 
sile strain on an ordered F-M can stabilize the FC-I as 
the ground state. This is because tensile strain naturally 
provides both BW reduction and suppression of superex- 
change and therefore stabilizes the FC-I phase. In par- 
ticular, the FC-I phase may thus be stabilized better (by 
suppressing other competing phases) in Lao.sCao.sMnOs 
by applying tensile strain. 

Effect of strain on magnetotransport. — Usually the 
maximum CMR temperature achievable by BW tuned 
phase competition is the Tc of the unstrained mate- 
rial. Additionally, the bicritical nature of the phase dia- 
gram and the proximity to the metal-insulator boundary 
needed for CMR conspire to keep the Tc quite low [33]. 
We now show that because strain affects different intrin- 
sic energy scales differently, it not only tunes phase com- 
petition, but also allows optimization of the competition 
between CMR temperature and %MR. 

Fig. 3(a) shows the resistivity, p, as a function of tem- 
perature for various tensile strain values applied to the 
unstrained F-M phase (blue dot in Fig. 2(b)). We see 



that increasing tensile strain causes a rapid rise in the re- 
sistivity maximum, that occurs at T ~ Tc, accompanied 
with a reduction of both the Tc and the temperature at 
the resistivity maximum (Tcmr)- The reduction in Tc is 
due to the approach to the F-M/SCO-I boundary by in- 
creasing the tensile strain, as depicted in the inset in (a). 
The reason for the increase in the resistivity maximum is 
the strain-induced enhancement of metal-insulator coex- 
istence at Tc as illustrated in (b). The color maps here 
depict the volume fraction of the insulating regions (red 
patches) embedded in the otherwise conducting back- 
ground at Tc. These insulating regions grow in volume 
with increasing strain and have short range (jr, 7r) CO 
correlations; the same correlations that one finds in the 
competing FC-I phase. The thermal fluctuations at Tc 
are typically dominated by the nearest metastable mini- 
mum, in this case the FC-I phase. Further, since strain 
controls the proximity to the F-M/FC-I boundary, in- 
creasing tensile strain makes the FC-I state progressively 
approach the energy of the F-M ground state, favoring 
greater insulating regions with short range (jr, 7r) CO cor- 
relations. If we start with other initial (unstrained) start- 
ing points, then the tensile strain can result in a F-M 
to SCO-I phase transition. The qualitative behavior of 
the magnetotransport is the same as in the F-M to FC-I 
case shown here. Typical data for transport and mag- 
netic field response near the FM/SCO-I phase boundary 
is shown in the Supplementary Information Section IV. 

In Fig. 3(c) we show the %MR, defined as 100 x [p(0) - 
p(h)]/p(0) and calculated at T = Tcmr, as a function 
of magnetic field for two cases. One case corresponds 
to %MR close to the F-M/FC-I phase boundary with 
(circles) and without (diamonds) strain. The other case 
is for %MR close to the F-M/SCO-I boundary. 
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FIG. 3. (a) Resistivity as a function of temperature for dif- 
ferent values of tensile strain. The strain values are marked 
on the individual curves, which are color coded with the ar- 
rows in the inset (which is a schematic of Fig. 2(a) showing 
only the tensile strain part) . Note that while the temperature 
at which the peak occurs (which is close to Tc) decreases, 
the peak value of the resistivity increases, (b) Real space 
snapshots of the charge ordering at Tc; the volume fraction 
of checkerboard charge ordered regions (shown in red) grows 
systematically with tensile strain. The labels A, B and C on 
them and in (a) denote the strain, temperature values where 
these were calculated. The magnetoresistance (%MR) as a 
function of magnetic field for ~1% strain on a system close 
to the F-M/FC-I phase boundary and ^2% on a system close 
to the F-M/SCO-I boundary. The unstrained values for both 
parent parameter points are shown for comparison. 

The amount of increase in the resistivity maximum 
with strain and the %MR depends on the domains of 
metastability of the competing phases, the type of the in- 
sulator and the largest tensile strain that can be applied 
before the system becomes insulating. It is clear how- 
ever that regardless of the nature of the metal-insulator 
phase competition, applying tensile strain yields an enor- 
mous enhancement of %MR (circles and stars) over the 
unstrained values (diamonds and triangles). 

Finally we demonstrate that compressive strain can 
drive insulators across the metal-insulator transition into 
the CMR regime. Fig. 4(a) shows the temperature- 
dependent resistivity with increasing compressive strain 
on SCO-I. At 0.8% strain, the insulator-to-metal transi- 
tion is accompanied by a CE-to-ferromagnetic transition. 
Increasing the compressive strain further causes a mono- 
tonic increase in Tc (also seen in the inset in (a)). The 
peak in the resistivity with increasing strain is also sys- 
tematically shifted to higher temperatures. In (b) we plot 
%MR at the temperature of the resistivity maximum as 
a function of strain. We also show the corresponding 
CMR temperatures (T CM r)- We find that %MR is re- 
duced upon increasing strain, as expected, but there is in 
fact an optimal region in which Tcmr can be increased 
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FIG. 4. Starting from an insulating parameter point, it is 
possible to apply compressive strain to convert it to a F-M; in 
(a) we show the resistivity of the resulting F-M as a function 
of temperature for various values of compression (color coded 
with the double arrows in the inset, which is a schematic 
of Fig. 2(c)). There is a pronounced peak which systemati- 
cally shifts to higher temperatures with increasing compres- 
sion, (b) Temperature at resistance peak Tcmr (left axis), 
which increases with increasing compression, and %MR (right 
axis) as a function of strain, which decreases as compression 
is increased. The %MR does not decrease drastically at first, 
and there is a region where Tcmr can be enhanced without 
losing the large %MR. 

without substantially reducing %MR. We have checked 
that all the above results survive in the presence of ran- 
dom spatial disorder which mimics A-site disordering. 

Conclusion. — We should stress a major difference be- 
tween the strain engineering considered here and the con- 
ventional approach using isovalent substitution. While 
both can be used to tune the bandwidth of a given ma- 
terial, uniform strain will not introduce the short-range 
disorder that naturally results from substitution. As a re- 
sult, the ferromagnetic Tc and Tcmr of strained materi- 
als should be higher than in their substitution-engineered 
counterparts, especially in low-dimensional geometries 
(i.e., films). In turn, the larger intrinsic bandwidth leads 
to larger %MR at smaller magnetic fields, since the ex- 
ternal field required to align the core spins is reduced. 
Our model of strain can be easily extended to incorpo- 
rate strain relaxation for modeling thick films, or to study 
how strain-induced domains and disorder impact magne- 
totransport. 

Our results indicate some promising future directions. 
First, we have demonstrated that strain gives access to 
a large phase space of new and accessible states for a 
given unstrained material. Second, strain need not in- 
troduce disorder, in contrast to chemical substitution. 
Finally, strain directly impacts orbital occupancy in a 
tunable way, and opens new possibilities for "orbitron- 
ics," or orbital-state sensitive electronics. 

We gratefully acknowledge support from DOE DE- 
FG02-07ER46423 (A.M.), DOE BES de-sc0005035 
(M.R.), NSF DMR-0907275 (N.T.), and the Center 
for Emergent Materials, NSF MRSEC DMR-0820414 
(W.S.C. and P.W.). 
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In this supplementary material we present the following: 

1. Model: First we give the details of the Hamiltonian and the parameters used in the main paper in Section I (A) 
and (B) respectively. We then discuss the scaling of the hopping and the superexchange parameters under strain in 
(C). There we argue why the change in the Mn-0 bond length change plays a dominant role in determining these 
scalings. 

2. Method: We present our method of solution in Section II. The definitions of the spin, charge and orbital structure 
factors and a brief account of transport calculations are discussed in Section III. 

3. Supporting data: In Section IV we show the magetotransport data (under strain) for a F-M parameter point 
close to the F-M/SCO-I phase boundary. When contrasted with Fig. 3(a) in the main text, they show that our 
conclusions regarding the strain response of magnetotransport are independent of the nature of the insulator-metal 
boundary near which the response is calculated. 

I. HAMILTONIAN AND PARAMETERS 
A. The manganite Hamiltonian 

Physical origin of the microscopic interactions: 




FIG. 1. Schematics of the relevant interactions incorporated in the manganite Hamiltonian: (a) The manganese level diagram. 
The largest scale is the Hund's rule scale that splits the up and down manifolds of the Mn 3d levels by Jh- Next in the 
energy hierarchy is the cubic crystal field, that splits that e g and the ti g states by /\cr ys {— lODq). For both the Mn +3 and 
Mn +4 states, three electrons occupy the ti g down manifold and are well-localized. The fourth electron, for Mn +3 , occupies the 
e g manifold, (b) As a result of placing a single electron in the degenerate e g manifold, the system spontaneously undergoes 
a Jahn- Teller distortion. The electron orbital-pseudospin operator couples to the local distortions Q2, Q3. (c) The effective 
antiferromagnetic coupling between two nearest neighbor Mn +4 atoms. The two p electrons on the oxygen connecting the two 
Mn sites can gain maximum kinetic energy by hopping to the two Mn sites if their ti g spins are antiparallel. This holds true 
even when the two Mn atoms are in the +3 state, except that the magnitude of this antiferromagnetic superexchange is lowered. 
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FIG. 2. (a) Experimental temperature (T) - phase diagram for ordered, half-doped manganites, from D. Akahoshi et a/. 1 , 
(b) Previous theoretical work, from K. Pradhan et al 2 : T/t — X/t phase diagram at x = 0.5 for ordered half doped manganites. 
There is a phase transition between the large X/t or small bandwidth (BW) SCO-I and the small X/t or large BW, F-M on 
increasing the ta- In the experiments, small ta implies small BW, which corresponds large X/t in theory. 



As shown in Fig. 1(a), the strong Hund's coupling and the crystal field ensure two things: first, that three electrons 
on the Mn site are well localized in the t^ g manifold, forming a S = 3/2 "core" spin that may be treated as a 
classical object, and second, that any electron occupying the e g level must have its spin parallel to the on-site tig 
core spin. The electron derealization happens through the e g manifold in the metallic state which is modeled by a 
two-band kinetic energy term with the on-site constraint that the itinerant spin orientation is projected on to the local 
"core" spin direction. This is the double exchange interaction. The spontaneous structural distortions arising from e g 
degeneracy-lifting when the e g manifold is singly occupied, as described in Fig. 1(b), is incorporated by coupling the 
e g electron's orbital pseudospin operator to the Jahn- Teller modes Q2 and Q3. Finally, the Mn-Mn superexchange, 
due to virtual exchange of the localized oxygen 2p electrons with the Mn e g states, is shown in Fig. 1(c). 

The Hamiltonian - Based on the preceding considerations, in this work we consider the microscopic two-band model 
for e g electrons with a strong Hund's coupling to tig core spins in a two-dimensional square lattice. The itinerant 
electrons are coupled to Jahn- Teller (JT) phonons and the core spins have a nearest-neighbor antiferromagnetic 
superexchange coupling between them. Explicitly, the Hamiltonian is given by 

H = E %p4 a *CjP*-JH^S i -<T i + J A FY l S i -S j 

(ij}a,a(3 i (ij) 

-A 2 {Q i ■ Ti + Q liPi ) + (jQl + PQl) 

+ ^2e a n itCt - fiN - h^2S zi (1) 

i,a> i 

where c and are the annihilation and creation operators for e g electrons and a, f3 index the two Mn e g orbit als 
(d x 2_ y 2 and <i 3z 2_ r 2, labeled a and b throughout). In the unstrained case, t^p are hopping amplitudes between 
nearest-neighbor sites with the symmetry-dictated form 3 : 

f/3, 

-t/y/3, 
t/y/3 

where x and y refer to the spatial orientation as in Fig. 1 in the main article. The e g electron spin operator is given 



4-x _ ,y = 
l bb — l bb — 

j_X iX 

l ab — l ba = 
L ab — L ba — 
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by <rf = Y^aa' c \acFaa' c ioicr'i where the r are the Pauli matrices. This spin is coupled to the local (classical) t^g 
spin Si via the Hund's coupling J#, and we assume Ju/t ^> 1. Finally, A is the coupling between the JT distortion 

Qi = (Q2i,Q3i) and the orbital pseudospin operator r- 1 — YTf c \ao^af3 c i(3°' This same X a lso sets the scale for the 
coupling between the breathing mode Qi and the on site charge density. K is the lattice stiffness. /?, the breathing 
mode stiffness, is assumed to be sufficiently large that this mode can be neglected. Finally, h is an external magnetic 
field, assumed to be in the z direction and coupled only to S{. We set t = 1, K = 1, and treat the Q i and the Si as 
classical variables which is a a good approximation for the manganites 3 . The chemical potential /i is adjusted so that 
the electron density remains at n = 1/2. 



B. Parameter values 

Fig. 2 shows the comparison of the experimental, temperature (T)-va phase diagram 1 and earlier theoretical work 
on the T/t -electron-phonon coupling (X/t) phase diagram employing the above model 2 . Both the experiment and 
theory are for ordered, half-doped manganites. The term 'ordered' implies alternating layers of A=Ba and A'=any 
Lanthanide (Ln), marked on the top of Fig. 2(a). Also the corresponding ta are shown on the lower axis. 

The value of ta controls the bandwidth (BW) in the experiments. The smaller the ta, the smaller the BW. In the 
theoretical work, shown in Fig. 2(b), the change of the BW is incorporated by the ratio X/t, where large X/t implies 
a smaller BW and vice versa. This not only serves as a benchmark of the model and the method of solution, to be 
described shortly, it also allows us to fix the Hamiltonian parameters by comparing the theoretical and experimental 
Tc values. 

For the electron-phonon and superexchange couplings, we choose the values of A and J to roughly reproduce the 
thermal magnetic ordering scales seen in experiments. Since we measure all parameters in units of the bandwidth 
(BW) £, the ratios X/t and J ft are what we fix. From Fig. 2, we find that for X/t ~ 1.5 and J/t = 0.1, T c ~ 0.045£. 
From this we obtain a Tc of about 225K, which is at least in line with the experimental value near the SCO-I/F- 
M phase boundary. In doing this conversion we have assumed t = 0.3eV. In general, there is a consensus that 
X/t ~ 1 — 1.6 and J/t ~ 0.05 — 0.1 is adequate for reasonably quantitative comparison with experiments for most 
doping values 3 and across many families of manganites. 



C. Scaling of hopping parameters and superexchange under strain 

(i) Strain puts constraints on the lattice parameter a which in turn changes both the Mn-0 bond length d and the 
Mn-O-Mn bond angle <\> as seen in Fig. 1 (b) in the main text. The p — d overlap integral V p da scales as d -3,5 , with 
d being the center of mass distance between the Mn and O atoms 4 . The Mn-O-Mn hopping is therefore proportional 
to Vp da / A and also depends on (j) as cos n ((/>) (as shown below). Here, A is the energy denominator that depends on 
the difference between the Mn and the O states involved. The dependence of the hopping integral on (j) is computed 
from the Slater-Koster tables 5 assuming the Mn-O-Mn bond to lie in the x-y plane. 

The exponent n is 3 for t aa , 2 for t^ and 1 for t aoi where a = d x i_ y i and b = 3 Z 2 _ r i. Thus the overall scaling for 
say, t aa would be d _7 cos 3 ((/>). Fig. 1 (b) in the main text, shows the relation between the lattice parameter (a), the 
Mn-0 bond length (d) and the Mn-O-Mn bond angle ((/>). Let us assume the strain imposes a change in the lattice 
parameter by Sa and this causes a change in d by Sd and <j) by S(j). Thus, t aa — » t aa = (d + Sd)~ 7 cos 3 (<j) + 5(j)). In bulk 
half doped manganites the value of (j) varies from 160° to 170° in SCO-I, Yo.sBao.sMnOs to F-M, Lao.sBao.sMnOs 6,7 . 
For these values of (j) and assuming small Sd/d and 50, we can linearize the expression for t aa in Sd/d and 5(j). This 
gives t aa ~ taa [l — 7^ — 3tan ((/>)£(/>] . The last term controls the extent to which the change in the unstrained Mn- 
O-Mn bond angle will effect the hybridization. For materials with <j) ~ 160°, the ratio |3tan(</>) • d^f^\ is 1 0.156 • djj\. 
To estimate this ratio we need typical data for perovskite structure under strain. Unfortunately, to the best of our 
knowledge, such data does not exist for the manganites. So as an estimate we assume a reasonable \5d/d\ ~ 0.01 at 
1% strain. Then it can be shown that for S(j) < 4°, the in-plane hopping increases (decreases) with Sd being negative 
(positive), while the bond angle dependence causes only quantitative change. Thus we neglect the S(j) dependence in 
our calculation for small strain. 

In Fig. 1 (b) in the main text, a = 2dsin(0/2). However in the spirit of above approximation we simply consider, 
a — 2d, or, 2Sd = 5a. So for the x-y plane, ey = Sd/d. Lao.7Sro.3Mn03 8 . Thus the hopping in the (x-y) plane scales 
with strain as t\\ ~ t (l — 7ey), while the out of plane hopping scales as t± ~ t (1 — 7e±). 

(ii) The scaling of the superexchange, J, with strain depends of the scaling of the fourth power of V v da under 
interatomic separation. Using considerations similar to above it is easy to show that, J\\ — )> J\\ ~ (Jy)(l — 14ey)) with 
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FIG. 3. The schematic of the SCO-I phase: (a) shows the planar checkerboard charge order (CO) with larger charge on the 
'yellow' sites and the CE type spin order (zig-zag ferromagnetic chains coupled antiferromagnetically) . (b) shows the alternating 
d x 2_ r 2 /d y 2_ r 2 orbital order (on the sites with larger charge density), (c) The color map of the magnetic structure factor, S(q) 
in the momentum space of our numerical data for the CE phase at low T. The q — (0,7r), (tt, 0), and (7r/2,7r/2) peaks are 
clearly seen. 



strain. Similar scaling has been observed on application of pressure in La2Cu04 9 . The anisotropy in the hopping is 
also reflected in anisotropic antiferromagnetic scales J\\ and J±. 



II. METHOD OF SOLUTION 



As mentioned above, we treat the core spins classically and the phonons in the adiabatic limit. Thus the solution 
of the Hamiltonian amounts to determining the configuration of the classical core spin and phonon background (for a 
fixed fermion density and temperature) that minimizes the free energy. To do this, we perform an exact diagonalization 
(ED) of the itinerant electron system for each configuration of the background spins and phonons. The Boltzmann 
weight of the combined electron, spin, and phonon system is calculated, and new spin/phonon configurations are 
sampled through classical Monte Carlo (MC). At a fixed temperature a Monte-Carlo system sweep consists of visiting 
every site of the system once in a sequential manner and performing the above mentioned update. 

The combined algorithm (ED+MC) is numerically rather costly, since the exact diagonalization must be performed 
at every step and the cost scales as 0(N 3 ) with N the number of lattice sites. Additionally, with a sequential system 
sweep, the cost of a Monte- Carlo system sweep scales as N 4 at each temperature. To reach reasonably large system 
sizes, we employ a recently developed variation 10 of real space ED+MC to that allows a linear scaling with the system 
size. This technique defines a region (cluster) around the site at which one attempts a MC update and accepts or 
rejects the update based on the energy change only in the cluster rather than diagonalizing the full system. For the 
calculation of observables, we diagonalize the full system, after the system has equilibrated using TCA. This adds 
only a few hundred full system diagonalizations to the computation cost. 

Within this 'traveling cluster approximation' (TCA), at each temperature, the computation cost for a system with 
N sites is O(Njj), where Nc is the fixed cluster size, and the cost of a full sweep of the lattice is linear in TV as 
opposed to A^ 4 . Using this technique we have accessed sizes up to 24 2 as opposed to the practical limit of ~ 8 2 within 
conventional ED+MC in 2 dimensions. In the present work on two dimensional systems, we employ a 8 2 moving 
cluster on system sizes between 16 2 (for most data presented here) to 24 2 . In studies such as this, where real space 
phase separation may play a vital role in determining transport responses, it is crucial to have access to system sizes 
that are large enough to capture coexisting phases on the lattice. 

We perform 4000 Monte-Carlo system sweeps at every temperature. Of these the first 2000 are used to equilibrate 
the spin and phonon variables. Every 100 th step from the remaining 2000 is used to calculate various quantities of 
interest. 

For calculating the response to an applied magnetic field, we repeat the above annealing process in the presence of 
the magnetic field. For example to calculate magnetoresistance at a given value of strain, we compare the resistivity 
calculated with and without a magnetic field. 
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FIG. 4. The resistivity vs temperature (a) and the response of the resistivity to magnetic fields (b). The data is shown for a 
F-M which is close to a tensile strain driven F-M/SCO-I transition. In the inset in (a), the location of Tco (yellow) and T c 
(rest) for the different values of strain, as indicated in main panel, are shown in the same color by the double headed arrows. 
We find that the strain and field response of the dc resistivity is qualitatively similar to the results for strain induced F-M/FC-I 
transition in the main text. Also shown in (b) for comparison, is the unstrained resistivity vs temperature at zero field (dashed 
line) . 



III. CHARACTERIZATION OF ORDERED PHASES & TRANSPORT 



A. Structure factors 



To characterize the phases as a function of strain, temperature, and magnetic field, we employ a number of static 
structure factors. The magnetic structure factor, defined by S(q) = ^2ij(Si ■ S j)e iq '^ Vi ~ Vj \ has a sharp peak at 
q = (0,0) for the ferromagnetic state and peaks at q = (0,7r), (7r,0), and (7r/2,7r/2) for the in-plane 'CE' type spin 
order. The long range charge order is characterized by the structure factor Dq(q) = ^^(Qi • Qj)e iq '^ Vi ~ Vj \ The 
in-plane checkerboard charge order is indicated by a peak at q = (ir, tt). The schematic of the SCO-I phase in shown on 
Fig. 3. (a) shows the a-b plane checkerboard charge order and the antiferromagnetically coupled zig-zag ferromagnetic 
chains, (b) shows the schematic of the alternating orbital order, (c) shows the color map of the magnetic structure 
factor, S(q) in the momentum space of our numerical data for the SCO-I order at low T. The q = (0,7r), (7r,0), and 
(7r/2,7r/2) peaks are clearly seen signifying the stabilization of the CE phase. 



B. Transport calculations 

The d.c conductivity o& c is estimated by the Kubo- Greenwood expression 11 for the optical conductivity. In a 
non-interacting system: 

*M = E(»« - np)^Ls(u - (ep - e a )) (2) 

a, (3 ^ 

The f a p are the matrix elements of the current operator, e.g., (ip a \j x \ilJp)i an d the current operator itself (in the 
tight-binding model) is given by j x = ita e Y^i,a( c l+a x,a c h^ ~ ^.c). The a are single-particle eigenstates, and € a are 
the corresponding eigenvalues. The n a = f(/i — e a ) are Fermi factors. 

We can compute the low-frequency average, <r av (/i, Acj, N) = (Acj) -1 J Q Aa; <r(/i, a;, N)duj, using periodic boundary 
condition in all directions. The averaging interval is reduced with increasing TV, with Au ~ B/N. Here the constant B 
is fixed by setting Acj = 0.06 at N = 40 2 . The d.c conductivity is finally obtained as cr^ c (/i) = limi y ^ 00 cr av (/i, B/L,L). 



6 



The chemical potential is set to target the required electron density n. This approach to d.c. transport calculations 
has been benchmarked in a previous work 12 . 



IV. TENSILE STRAIN INDUCED METAL TO INSULATOR TRANSITION ACROSS THE F-M/SCO-I 

BOUNDARY 

Here we show that qualitatively our conclusions are independent of the choice of the metal-insulator boundary near 
which we calculate magnetotransport. For this we present the magnetotransport for a F-M with (X/t = 1.5, J/t = 0.1) 
which is close to tensile strain driven F-M/SCO-I boundary. Fig. 4(a) shows that the enhancement in % MR with 
tensile strain is similar to that seen in Fig. 3(a) in the main text. Also shown in Fig. 4(b) is the magnetic field 
induced insulator to metal transition for the case of 2.2% tensile strain. Fig. 3(c) in the main text, are constucted 
from similar data. 
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